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Abstract 

Some speculative preliminary ideas relating matrix theory and cosmol- 
ogy are discussed. 



1 Slightly extended version of the talk given by G. W. Gibbons at the Mitchell Institute 
Conference on Cosmology 
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1 Motivation 

This is a report on some on-going work in which an attempt is made to ex- 
plore how to incorporate the basic ideas of cosmology into M-Theory. It may 
be seen either in the context of much recent work on time dependent back- 
grounds in String Theory, or in its own right, as a speculative approach to cos- 
mology aimed at ultimately taking us beyond the standard Friedman-Lemaitre 
paradigm. More concretely, our motivations are 

• The BFSS matrix model is claimed to provide a fundamental quantum 
mechanical description of "M-Theory" 2 . 

• It replaces conventional spacetime concepts, such as commuting coordinates, 
with inherently non-classical notions such as non-commuting coordinates. 

• It should, therefore, surely have something deep to say about the structure of 
the universe. 

2 For a useful review of M(atrix) theory see e.g. [3]. 



2 



• In particular one should be able to use it to address such issues as the existence 
and significance of such things as "The Wave function of the Universe" . 

In this talk we present some rudimentary and very preliminary ideas aimed at 
understanding how we should think about cosmology in the language of matrix 
theory. It is a report of work done partly in collaboration with S. Alexander. 
The only previous work on this topic known to us is that of Alvarez and Meessen 

Q. 

2 Matrices and DO-particles 

One may view the matrix model in two slightly different ways. 

• Following BFSS, as the limit N — > 00 of a super-quantum-mechanics of 9 
N x N Hermitian matrices. 

• Following earlier work by de Wit, Hoppe and Nicolai jU, as a regularization 
of the super- membrane of 11-dimensional supergravity. 

• Both approaches lead, because of the high amount of super-symmetry, to 10- 
dimcnsional super- Yang-Mills with gauge group G = U(N) and fermions in the 
adjoint representation reduced to one spacetime dimension. In Coulomb gauge, 
one replaces the u(N) valued connection one- forms (x, t) by their 9 spatial 
components Ai{t) which are the 9 Hermitian matrices X l (t) of the model. 

• From the membrane point of view one passes to light-cone gauge and the 
X 1 represent the 9 transverse components of the membrane coordinates. The 
residual bosonic gauge-invariance consists of sdiff(S2), area preserving diffeo- 
morphism of the membrane 2-manifold S2. The Lie algebra of sdiff(S2) is well 
known to coincide, in some sense at least, with limjv^oo u(N). 

• From the 10-dimensional String Theory point of view one should regard the 
X 1 as representing the 9 non-commuting position coordinates of N DO-branes, 
the locations of the ends of fundamental strings. 

• Clusters of large numbers of DO-particles are described by classical solutions 
of 10-dimensional Type IIA supergravity theory. 

The BPS states correspond to electrically charged singular extreme "black hole" 
hole solutions where the electric Ramond-Ramond charge couples to a gravi- 
photon field of the 10-dimensional Type IIA super-gravity theory which may be 
obtained by dimensional reduction of the 11-dimcnsional super-gravity theory. 

The classical solutions, which describe k separated clusters in force balance 
may be lifted to 11- dimensions where they have the structure of k singular 
vacuum pp-wave solutions moving parallel to each other. The pp-waves are 
sometimes described as 11-dimensional gravitons, but this is not really accu- 
rate, because even in 11-dimensions, the solutions have distributional sources. 
A better description is as lightlike cylinders extending along the 10'th spatial 
dimension. 
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3 Classical Matrix Theory 



One approach to Matrix Cosmology is via the classical equations of motion [3] . 
We shall briefly describe this, since it was the original approach that we adopted, 
but later we will, for reasons to be explained, abandon it for a rather different 
picture. 

The basic classical equations of matrix theory are 



(1) 



where X % are n N x N hermitian matrices, the index j is summed over and 
A = A^- is a possible cosmological or mass term [5] El El El- 

In the BFSS model n — 9, and one is looking at N DO-branes but these 
equations have been studied more widely as a reduction of U(N), or, if they 
are taken traceless, SU(N), Yang- Mills theory to one time dimension. As such, 
there is some evidence for chaotic behaviour. 

Note that 

• If A = 0, we have a Galilei invariant system 

• If A ^ 0, we have invariance under one of the two Newton-Hooke groups, of 
the two non-relativistic contractions of the De-Sitter (A > 0) or Anti-de-Sitter 
(A < 0) groups. A description of these groups and their transformation rules 
together with an account of their significance for Newtonian cosmology with a 
cosmological constant are given in [7j. 

• If one thinks in terms of a mass term, then a positive mass squared corresponds 
to negative cosmological constant and a tachyonic mass term to a positive cos- 
mological constant. 

• In the case of the BFSS model, the equations of motion must be supplemented 
by a constraint on the initial conditions which arises from the Gauss constraint 
of the gauge theory 

[X\X j } = 0. (2) 



d 2 X l 



[X j , [X^X 1 ]] = \x l 



4 Newtonian Cosmology 

Since our equations lack manifest covariance, the most helpful analogy is with 
elementary Newtonian Cosmology in 3 spatial dimensions The brief presentation 
of Newtonian cosmology which follows may be unfamiliar, but it is completely 
equivalent to more conventional accounts in the literature. The generalization 
to other space dimensions, and indeed to other force laws, is trivial. 
Newton's equations of motion for k gravitating particles are 



y m a mb{vb - r„) , , 
G — : - h Am a r Q . 



(3) 
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• As with the matrix model, we either have Galileo invariance (A = 0) if the 
cosmological constant vanishes, or if it does not we have Newton-Hooke invari- 
ance. 

• In order to incorporate the Cosmological Principle we make a nomothetic 
Ansatz 

r«(t) = a(t)Xo, (4) 

where the so-called co-moving coordinates x a are independent of time. 
The homothetic ansatz leads to two conditions. 



3 + A. (5) 



• Raychaudhuri's Equation 

a, p 
a 3a 

This is the usual equation of motion for the scale factor of an expanding universe 
with a cosmological term and pressure free fluid. In what follows we set A = 
for simplicity. It is easy to adapt the discussion to the case A ^ (see [7]). 

• The co-moving coordinates must constitute a Central Configuration, i.e. a 
solution of 

-m a x a + ^G |xg _ x&|3 =0. (6) 

The quantity /x is a constant. Central configurations are extrema of an auxiliary 
potential 

Recently, with Battye and Sutcliffe |S] , one of us has carried out an extensive 
numerical investigation of central configurations which are minima or ground- 
states of this potential for up to 10 particles. In the case of equal masses 
m a = m, Va, the conclusion is that the minima correspond to a spherical ball 
of particles of uniform density j^q- In other words, if N is the total number of 
particles and the radius p is defined by 

GNm _ pp 

then one finds a uniform density of particles inside the radius p and almost no 
particles outside that radius. The interpretation of J3J should be clear. It is 
well known that a spherical shell of matter exerts no force on particles inside 
it but an attraction on particles outside given by the total mass of the shell. 
The left hand side of JHJ is the Newtonian attraction due to the total mass Nm 
interior to radius p on the thin shell of particles at p. The right hand side is 
the repulsive pseudo-cosmological force on the thin shell which is proportional 
to the distance p. We use the term 'pseudo-cosmological' to alert the reader to 
the fact that we get such a term even if the cosmological constant A = 0. It 
really arises from the inertial term in the Newtonian equation of motion ©. 
Note that the proper radius of our ball is time-dependent and given 

\v\=R = a{t)p. (9) 
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A ball of uniform density is exactly what one expects on the basis of the 
usual pressure free fluid model. Thus this, slightly unconventional, approach to 
Newtonian Cosmology reproduces all of the standard features without making 
arbitrary assumptions about fluids, rather these assumptions are derived from 
the model. 



4.1 Quantum Newtonian Cosmology 

In order to prepare ourselves for Quantum Matrix cosmology, it may be worth 
pausing to recall that one can obviously construct a Wave Function for Newto- 
nian Cosmology in the framework of non-relativistic quantum mechanics. This 
may not often be done in discussions of Quantum Cosmology but it is entirely 
straightforward and elementary. All that one needs is a solution of the multi- 
particle Schrodinger equation 



with * = tf(r«), V 2 a = -jL and 



dt ^— ' 2m a 

m a m b 



4.2 WKB Approximation 

We consider a potential V — V(r a ) which is homogeneous of degree n. For the 
case of Newtonian gravity without a cosmological term corresponds to n = — 1. 
The equation of motion is 

dV 

m a v a = (12) 
or a 

The homothetic ansatz is 

r a = a(t)x a , (13) 
where the co-moving coordinates x a constitute a central configuration satisfying 

^m a x Q = +— 14) 
3 <9x a 



and the scale factor satisfies the Raychaudhuri type equation 

a 1 -»a = -^. 

3 

with first integral or Friedmann equation 



(15) 



la 2 + -^a n =-, (16) 
2 n 3 2' V ' 

where k is a constant. Taking the dot product of x a with (|14f> and using Euler's 
theorem gives the Virial Theorem 

|£m a x2=nV(x ). (17) 



G 



The conserved energy is 



2 

At the JWKB level, the wave function is 



H=|$> a x2. (18) 



Wwe", (19) 
where S is the relevant solution of the Hamilton- Jacobi equation 



1 fdS\ 2 rr/ . dS 



2m a \dr a J ' dt 

In our case the relevant solution is 

v a m a 2 \ -> . m a 2 .if 

4.3 Hartree-Fock approximation 

Here we suppose all masses equal m a = m and replace the full wave function 
&(r a ) by the product 

*(r a )<xJJ*'(r a ), (22) 

a 

where \I/'(r) satisfies 

*^ = -^V 2 *' + mi7*', (23) 
at 2m 

2 

a? 7 

V 2 C7 = 4irGm\^'\ 2 . (24) 



with W = *'(r,t), V 2 = |^ and 



The time-independent Schrodinger equation coupled to Poisson's equation has 
been studied in a different context where it is referred to as the Schrodinger- 
Newton equation [§]. 

In our case, we assume that 

U = F(t)r 2 , = A(t)e% , (25) 

We find from the Poisson's equation that 

One readily sees that one must have 

S = B{t)r 2 , (27) 
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with 



SB _ A 
m A 



(28) 



Moreover if A = a ? , then the scale factor a(t) satisfies the Raychaudhuri 
equation 

4?rGm (29) 



The action is given 

and the wave function by 



3a 2 ' 

a 1 " 2 

S = -m-r 2 , 
2 a 



0,2 



(30) 



(31) 



In this case the Hartree-Fock approximation gives a version of the WKB wave 
function corrected by the prefactor —3-. 



4.4 Normalization and Energy 

In order to model the central configurations described earlier which have a 
finite number of particles, we need to use a normalizable wave function. We 
use the Hartree-Fock wave function ^'(r) of Q31(l in the region r < R and take 
^'(r) = for |r| > R. This is an exact solution of the Schrodinger equation 
(|23|l within each region, but fails at the surface |r| = R. We ignore this issue 
here. 

The normalization integral of ^'(r)! 2 has support in a ball of proper radius 
|r| = R. The norm must be time-independent, so we need to take time dependent 
R(t). In fact we need to take 

R(t) oc a(t) (32) 

This is consistent with the classical analysis giving Hubble's law ©. 

The classical problem of Newtonian cosmology has a conserved energy, and 
we should check the energetics of our quantum-mechanical model. The energy 
of a single particle wave function is 



d A r 



2m 1 1 



2l,Tr'|2 



mF(t)r 2 \^ 



oc 



-B(t) 2 + F(t) 



■ R(t) 

A(t) 2 I d 3 rr 2 



1 . ,0 2nGm 
The Raychaudhuri equation (|29() has the first integral 



-At)' 



inGrri 1 
lo(t)" = 2 ' 



(33) 
(34) 
(35) 

(36) 
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so we have 

,_, 1 , 2nGm ,„„. 

This is a constant, as desired, plus a i-dependent error term which we attribute 
to the sharp cutoff in the wave function. 

Since the quantum-mechanical energy l|33|) is strictly positive, we must choose 
the k > solution of (|29|l with large t behavior a(t) ~ kt. The error term above 
vanishes at large t. The quantum mechanical model is thus consistent for an 
"open universe" . 



4.5 Wick Rotated Newtonian Wave function of the Uni- 
verse 

We have constructed an approximate wave function of a simple Newtonian uni- 
verse whose WKB approximation gives a classical solution of Newton's equations 
of motion representing an expanding gas of point particles. The aim of quantum 
cosmology is to derive this wave function, and hence the initial conditions for 
the universe from some more fundamental assumption, such as the No-Boundary 
Proposal of Hartle and Hawking. We shall not dwell on this in detail here but 
content ourselves with the following, possibly suggestive, remark. If we take the 
simplest (Einstein-de-Sitter) solution for the scale factor a(t) oc is, we have 

f'oce^ 12 ^. (38) 

Curiously, this Euclidean wave function, strictly speaking a solution of the dif- 
fusion equation rather than the Schrocdinger, will be normalizable with respect 
to integrations over the positions if we Wick rotate, i.e. set 

t = -it, (39) 

with the imaginary time coordinate r being real and positive. One might spec- 
ulate that this normalizability of the Wick-rotated Newtonian wave function is 
related to Hartle and Hawking's path integral approach to the wave function of 
the universe. 



5 Homothetic Matrix Cosmology 

After preparing ourselves with a brief excursion into Newtonian Cosmology, we 
return to the matter at hand. In matrix cosmology is it also natural to begin 
by making a homothetic ansatz 



where M l are independent of t. 
Substitution leads to 



X l = a{t)M\ (40) 



-3 - s = M, 41 
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fiM i + [NP , [M j , M 1 ]] = 0. (42) 

The idea is to interpret (|41|l as the analogue of Raychaudhuri's equation in 
cosmology and (|42|) as the analogue of the equation governing central config- 
urations in Newtonian cosmology 8 or monopole scattering Essentially 
the same equation arises in supersymmetric J\f = 1* gauge theories when one is 
looking for vacua or ground states ^21- For that reason we shall sometimes refer 
to solutions of (|42|l as vacua. Similar equation also arises in certain solutions 
describing spinning membranes |13| . 

We begin by looking at Raychaudhuri's equation l|41|) . It has a first integral 

a\ k iia 

^ + ^V +A > (43) 

where A; is a constant of integration. This is the analogue of the Friedman 
equation in standard cosmology . As well as a standard cosmological constant 
or dark energy contribution given by A we have an exotic matrix contribution 
to the energy density given by 

PM= 1^G- (44) 

If fj, is positive this energy density increases as the universe expands, indicating 
that the pressure Pm has the opposite sign to the energy density. In d spatial 
dimensions this would lead to 

d + 2 fAKS 

Pm = —Pm ■ (45) 

From the point of view of matrix theory the most natural choice for d would 
be 9. Arguing by analogy with the Newtonian case one might then regard a(t) 
as the scale factor in Einstein conformal gauge. This then leads to 

Pm = ~Pm- (46) 

Later we shall compare this with a model based on a supergravity solution 
representing a gas of expanding DO-branes. 

As an example of the general theory, consider 3 NxN matrices M % , i — 1,2,3 
providing an N dimensional representation of su(2), 

[M 1 ,M 2 ]=iM 3 etc. (47) 

This solution actually describes an expanding spherical membrane (see |2] and 
references therein). From (|42|) we find 

pi = -2, (48) 

which implies a negative energy density and positive pressure. This looks rather 
unphysical and so we turn to an anisotropic model. Recall that 

9 = 3 + 3 + 3, (49) 
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and take 3 mutually commuting sets of such matrices, each with its own scale 
factor a(t), b(t), c(t) say. If M 1 , M 2 , M 3 are taken to be diagonal and A > we 
shall get exponential expansion for the scale factor a since 

a = \a . (50) 

On the other hand, we can have X 4 , X 5 , . . . X 9 oscillating. In other words 



3 directions exponentially expand 



(51) 
(52) 

This phenomenon is closely related to the well-known chaotic behaviour of 
Yang-Mill reduced to one time dimension and zero space dimensions. 



6 directions oscillate 



5.1 Chaos 

In standard su(2) Yang-Mills, one may assume that the connection is 

A = X i dx i = uit^dx + v{t)r 2 dy + w(t)T 3 dz , (53) 

where t 1 are Pauli matrices. The equations of motion derive from the La- 
grangian 

L = i(w 2 + v 2 + w 2 ) - V(u,v, w), (54) 

with 

V(u,v,w) = -(u 2 v 2 +v 2 w 2 +w 2 u 2 ) . (55) 

The non-negative potential V has three commutative valleys along the three 
orthogonal coordinate axes in u,v,w space for which V vanishes. Studies of 
the motion ^3] show that the representative particle rattles along each valley, 
eventually returning to the origin and rattling along another valley. This is 
rather reminiscent of the behaviour of the three scale factors a, 6, c of a Bianchi 
IX chaotic cosmology of the type originally studied by Misner ^H] . Introduction 
of a positive cosmological term, (i.e. a tachyonic Higgs mass) leads to the 
eventual escape of the particle in one direction, provided A exceeds a certain 
threshold. 

An obvious extension of this idea is to consider an ST/ (2) x SU(2) x SU(2) 
model in which just three directions expand exponentially and the other six 
remain bounded. 



6 DO-particle Cosmology 



According to Type IIA ten-dimensional Supergravity, DO-branes correspond to 
extreme black holes, a static configuration of k clusters depending on a harmonic 
function H on E 9 

k 



H = 1 
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where p, a is quantized, being proportional to N a for N a DO-branes located at 
positions x a . 

The moduli space is clearly given by k points in E 9 . The slow motion is 
governed by a metric induced from the De-Witt metric of the Type IIA action. 
Long ago, Shiraishi showed that this metric is flat [T2|. In other words, there are 
no velocity dependent forces quadratic in velocity. Thus one might anticipate 
that a cosmology of DO-branes should expand freely like a ten-dimensional Milne 
model which has 

a(t)<xt, p + 9P = 0. (57) 

This is contradicted by some exact solutions of Type IIA found by Maki and 
Shiraishi |17j some time ago, following earlier work by Kastor and Traschen. In 
these one has 

a(t)octi, p = P, (58) 

which corresponds to 'stiff matter', for which the sound speed is that of light. 
Neither of these two equations of state coincides with that given by the homo- 
thetic matrix model of the previous section. One might try, as was suggested 
to us by Justin Khoury to 'save appearances' by passing to string conformal 
gauge, which might seem more appropriate for the homothctic matrix model 
scale factor. However, in this gauge, we would have 

P = -\p, (59) 

which does not coincide with the P = —^p which we obtained from the homo- 
thetic matrix model. In hindsight this is perhaps not so surprising, since the 
matrix configuration leading to this equation of state were rather delocalized 
and in fact corresponded to extended objects. For well localized D-particles with 
mutually commuting coordinates the quartic term in the matrix Lagrangian is 
effectively zero and therefore the dominant interaction between D-particles will 
come from a 1-loop correction which generates a term v 4 /r 7 in the 2-body 
Lagrangian. It would be interesting to see whether Newtonian approach to 
cosmology based on those interactions could generate Maki-Shiraishi solutions. 
We shall not explore that in detail here but rather we shall describe the Maki- 
Shiraishi metrics and their 'hidden supersymmetry'. There is no denying that a 
cosmology made up entirely of DO-branes, with no anti-DO-branes might not be 
thought of as being very realistic. Nevertheless, the solutions we are about to 
discuss do exhibit some extremely interesting features which we hope contains 
lessons for future, more realistic models. 

6.1 Maki-Shiraishi metrics 

Maki and Shiraishi |17| considered as a Lagrangian in n+1 spacetime dimensions 
for gravity plus a two-form plus a scalar 

L = R — (V^) 2 -e^V-(n-l)e^*A. (60) 
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where a and b (and A) are constants. 
They sought solutions of the form 



ds 2 = -H ^T^dt 2 + a 2 (t)H^T^dx 2 . (61) 

H = 1+ a(tT- 2 +i ^(n- 2)|x - x a |«-2 ' (62) 
e& + = a(t) p H'^T^ , (63) 



n-1 dt /, 1\ 
A= V 2(n-2 + oa) ^)T( 1 -Hj' (64) 

with F = cJA. Maki and Shiraishi found various solutions. The time dependence 
of the scale factor depends on the particular solution, for us the relevant one 
satisfies n = a 2 = 9, p = 2a 2 , A = and 

a(*)=(£f. (65) 
If no DO-branes are present the background metric is 

ds\ Q = -dt 2 + a 2 (t)efoc 2 , (66) 

with 

g = e <t> = t%. (67) 

This is just what one expects for gravity coupled to a massless scalar field 
which behaves just like stiff matter. From the point of view of string theory, we 
see that we have a time dependent string coupling constant g which increases 
with time from a zero value at the Big Bang. In other words 



late times < — > strong coupling (68) 



early times < — ► weak coupling (69) 



This feature remains true if DO-branes are present. First note that one may 
take a(t) to be constant by making the change 

t->t + t a , (70) 

and letting to — > oo. One then obtains the static multi-brane solutions 

ds 2 1Q = -H-idt 2 + hUx 2 , (71) 

with 

H - 1+ t^r ,72 > 
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g = e* = H*, (73) 

and 

A = i 1 i) dt (74) 

Now let us restore the time dependence. One finds, setting to = 1, 

ds 2 w = -H-idt 2 +t%H^dy? , (75) 

with 

H = 1 + y t^r L — - r (76) 

^7t-|x-x a | 
g = = tiffi , (77) 

and 

Evidently the general time-dependent solution represents a gas of DO-branes in 
10 dimensions in a background time dependent dilaton field. Note that 

• the non-interacting gas of DO-branes does not affect the law of expansion. 

• while the solution at large distances is time-dependent, with the physical 
separation of the DO-branes increasing with time, near each singularity, i.e. as 
x — > x a , the solution is effectively static. 

7 Lift to 11 dimensions 

The Maki-Shiraishi metrics are clearly not supersymmetric, i.e. they are not 
BPS because they are time dependent. Every Lorentzian spacetime admitting a 
Killing spinor field must admit an everywhere non-spacelike Killing vector field. 
However if one lifts the solution to eleven dimensions using the uplifting formula 

dsl-t = e~* dsl + e~ (dz + 2Af , (79) 
where z is the eleventh coordinate, something interesting happens. 

7.1 The background 

Using the uplifting formula one finds that the background <|66II67|I becomes 

ds\ x =t~%(-dt 2 + t§dx 2 ) + t^dz 2 . (80) 

If one defines 

T=lti, (81) 
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we get 



dsn 



-dT l 



T 2 dz 2 + dx 2 . 



(82) 



This is flat space E 9 xE 1,1 in Milne coordinates. Of course if the tenth coordinate 
z is taken to be periodic, as it would be on the M-Theory circle, then we shall 
get the usual orbifold singularities and non-Haussdorf behaviour associated with 
Misner spacetime |18) . 



7.2 Lifting the general solution 

We define 

dt = tidt 



=> t = —t~ 
16 



We obtain 



Now let 



One gets 



ds^ = dx 2 + t~H(dz + 2 A) 2 



dz = dz + — . 



ds^ = dx 2 + t~H(dz) 2 - 2dzti 



dt 2 
tfH 



(83) 
(84) 

(85) 
(86) 



This looks complicated, but if we define a time independent harmonic func- 



tion 



and set 



H 



Mo 



x - x a 



T = -a 8 
9 ' 



t + t 



x° = Tcosh 



9t 



.10 



= T sinh 



8z 

9h 



x° ± x 10 



with 



we have 



dt 8 dx + 



dsn = dx 2 — dx + dx + H 



dx^ 



(87) 

(88) 

(89) 
(90) 

(91) 
(92) 



Note that 

•This is a pp-wave whose profile depends on light-cone time x + . 
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• The solution is nevertheless boost-invariant, the scalings 

x + -> \x + , x" -> A _1 x _ (93) 
with A 6 R\0, leave the metric invariant. 

• Reduction on the boost Killing vector gives the 10-dimensional solution. 

• The 11-dimensional solution is BPS, it admits a covariantly constant Killing 
spinor but this is not invariant under boost and hence the 10-dimcnsional solu- 
tion is not BPS. 

8 Conclusions 

• Homothetic solutions of classical matrix theory resemble expanding universes 
but do not really capture the cosmology of DO-branes. 

• Exact supergravity 10-dimensional Type II A solutions for expanding universes 
of DO-particles are available. 

• Lifted to 11-dimensions they are vacuum pp-wave solutions with time depen- 
dent profile and hence BPS. Their reduction to ID-dimensions is on a boost 
Killing field and hence they are time dependent and non-BPS in 10 dimensions. 

• It seems that Quantum Mechanical matrix theory in a suitable limit captures 
the behaviour of the classical super-gravity solutions. 

• The status of the 'Wave function of the Universe' remains unclear. 
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